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The three-gluon and ghost-gluon vertices of Landau gauge Yang-Mills theory are investigated in 
the low momentum regime. Due to ghost dominance in the infrared we can use the known power 
law behavior for the propagators to determine analytically the complete momentum dependence 
of the dressing functions. Besides a uniform, i. e. all momenta going to zero, divergence, we 
find additional singularities, if one momentum alone goes to zero, while the other two remain 
constant. At these asymmetric points we can extract additional infrared exponents, which cor- 
roborate previous results and expand the known fixed point solution of Landau gauge Yang-Mills 
theory, where the uniform infrared exponents for all vertex functions are known. Calculations in 
two and three dimensions yield qualitatively similar results. 

We find several dressing functions diverging like {p^Y^^^, if only the gluon momentum p goes 
to zero. Of these many are longitudinal and do not contribute to Dyson-Schwinger equations. The 
divergent transversal parts are additionally suppressed by the corresponding tensor. The longitu- 
dinal dressing function of the ghost-gluon vertex behaves similarly, when the gluon momentum 
becomes small compared to the ghost momentum, whereas all its other dressing functions van- 
ish at the asymmetric points. The uniform momentum dependence of the three-gluon vertex is 
determined as {p"^)^^^, while the ghost-gluon vertex stays finite in this limit. 
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1. Infrared Behavior 

Many non-perturbative aspects of Yang-Mills theory are encoded in the infrared (IR) behav- 
ior of its Green functions. For very small momenta below the intrinsic scale Kqcd they can be 
described by power laws. This scaling solution [1] is in accordance with the scenarios for confine- 
ment of Gribov-Zwanziger and Kugo-Ojima. The general expression for the so-called IR exponent 
of an arbitrary Green function with m gluon and 2n ghost legs under the assumption that all mo- 
menta go to zero uniformly is [2] 

&2n,m = {n-m)K+{l-n)(^-l}j. (1.1) 

This formula is valid in J = 2, 3,4 dimensions for the corresponding values of K. As it is clear that 
we can have more than one independent momentum for vertex functions, the question arises what 
happens when only one of these goes to zero? Interestingly it turns out that additional divergences 
can occur that do not change the uniform solution [3]. 

For the calculation of the three-point vertices we used the known power laws for the propaga- 
tors and the following truncations for the Dyson-Schwinger equations (DSEs), which is motivated 
by ghost dominance in the IR [1, 4]: 




For the dressed ghost-gluon vertex we used the bare one, justified by a simple argument of 
Taylor [5], which is supported by lattice simulations [6] and the DSE solution [7]. The truncation 
considers only the first order of the skeleton expansion. For more details see [8]. 



2. Three-point integral 

The ghost-triangle diagram of the three-gluon vertex is decomposed in the IR into ten ten- 
sors T'^vpipijPijPj) and ten corresponding dressing functions Ei{pi,p2,P3). The applied tensor 
decomposition [9] reveals that only ten instead of the expected fourteen scalar functions and 
tensors z'^yp are necessary. The former consist of massless three-point integrals, 

^ ' (2 1) 

{2nr{{q+p,ry'{{q-P2rr-{^^r'' ' 

where Vi , V2 and V3 are non-integer numbers. We employed the Negative Dimensions Integration 
Method (NDIM) which yields a full analytic solution in terms of Appell's series F4. Using several 
different analytic continuations we can calculate the momentum dependence of the diagram for 
ai^bitrai^y d. As the variables of the Appell's series are the momentum ratios p\/ p^ and p\/ p\ we 
plot the scalars Ei from the tensor decomposition as functions of p\ and P2 with p\ fixed. The 
plotted region in fig. 1 represents the accessible ratios for Euclidean momenta and is restricted by 
momentum conservation. The ghost-gluon vertex can be decomposed similarly into two parts. 
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Figure 1: Left: A double-logarithmic plot of the scalar functions £, of the three-gluon vertex when p\ 
becomes small and p\ = P2 is held fixed. Three cases can be distinguished: Five scalar functions diverge 
like {p\y^^^, three stay constant, two vanish like {piY'l^^'^^ and (pi)^^^'^ respectively. Middle and right: 
The scalar functions £5 and £10 with p\ = I. The kinematic singularities at the asymmetric points, where 
p\ or p\ go to zero, clearly dominate the structure. 

3. Three-point vertex functions 

As an example for the kinematic dependence of the dressing functions we show £5 and £"10 
in fig. 1. The behavior when one of the three momenta gets small compared to the others, i.e. 
the region around the asymmetric points (0, 1), (1,0) and (00,0°), is of special interest. This case 
corresponds to the emission/absorption of a soft gluon. We extracted the exponents in d dimensions 
and indeed found kinematic singularities, which are in agreement with the power counting analysis 
in ref. [3]. In fig. 1 the results in four dimensions are given. 

To verify the self-consistency of the assumption of a bare ghost-gluon vertex we calculated 
its momentum dependence for uniform scaling and the case of only one small momentum. The 
results for the former clearly support the bare version, whereas those for the latter show that the 
structure of the ghost-gluon vertex is richer than expected and features a divergence (1 — 2 K") for the 
longitudinal scalar function when the gluon momentum vanishes. However, this is no contradiction 
to Taylor's ai^gument because this scalar function is multiplied with the soft momentum and the 
divergence is suppressed. 
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